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ABSTRACT 

We present deep 850^m imaging of the z=0.773 strong lensing galaxy cluster RCSJ022434-0002.5 
from the Red-Sequence Cluster Survey (RCS). These data are part of a larger submillimeter survey of 
RCS clusters, with SCUBA on the JCMT. We find five objects at 850/Ltm, all of which are also detected 
at either 1.4-GHz, 450/im or both. The number density of objects in this field is in general agreement 
with the blank-field source counts; however, when combined with other cluster surveys a general 
tendency of cluster fields towards higher submm number densities is seen, which may be the result 
of unrecognized submillimeter luminous cluster galaxies. Primarily employing optical photometric 
redshifts we show that two of the five submillimeter galaxies in this field are consistent with being 
cluster members, while two are more likely background systems. 

Subject headings: (ISM:) dust, extinction - galaxies: clusters: individual (RCSJ022434-0002. 5) - galax- 
ies: evolution - galaxies:starburst - infrared: galaxies - radio continuum: galaxies 



1. INTRODUCTION 

The study of galaxy clusters encompasses almost ev- 
ery aspect of cosmology, from the determination of fun- 
damental constants to the influence of the local envi- 
ronment on galaxy evolution. Yet a solid understanding 
of the process of cluster evolution has remained elusive. 
In a hierarchical universe clusters begin to form at very 
high rcdshift, then grow and evolve through cluster merg- 
ing and accretion of groups. Understanding this pro- 
cess,which may continue to significantly lower redshifts 
than the original epoch of formation, and its importance 
in driving cluster galaxy evolution, will be a major step 
towards clarifying the complex physical details of hier- 
archical structure formation, at the level of individual 
objects. 

By a redshift of z ~ 1 galaxy clusters show al- 
ready significant differences from their local descen- 
dants , indicative of younger and more active sy stems 
(e.g..lButcher fc Oemlerlll984Hvan Dokkum et al'] | 199fi 



| 200ii iDwa rakana th fc Owed 119991 ISmail et all ll99S ; 
IBest et al .1120021) . As in the field, m uch of this activ- 
ity is dust enshrouded l)Poggianti et alll999HSmail et alJ 
1999) and conclusions drawn solely from optical studies 
are potentially misleading or at best incomplete. Indeed, 
recent observations of moderate redshift clusters (z ~ 
0.2-0.4) with ISO at 15/im, which directly samples hot 
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dust, have detected dusty star bu rts l)Fadda et alJ 1200(1 
IDuc et al.H2002HCoia et al.ll2004[) . and it is estimated for 
one cluster that perhaps 90% of the star forma tion is dust 
enshr ouded and undetected in optical studies ijDuc et alJ 
2002). The ISO results are parti cularly interesting since, 
at least for one cluster, A1689 (|Fadda et al.ll2000lh the 
fraction of dusty galaxies is greater than in the field, im- 
plying the cluster environment is somehow responsible 
for triggering the intense star formation. 

Thus far, extreme star-formation activity in clusters 
(such as that seen in ultraluminous infrared galaxies 
[ULIRGs] at high and low redshift) has not been unam- 
biguously detected; even the systems seen by ISO are lim- 
ited to rates of a few tens of solar masses per year (assum- 
ing no contribution from active galactic nuclei [AGN]). 
The exception to this is the handful of central cD galax- 
ies detected at 850M m primarily during surveys for faint 
background sourc es (Ed ge et alJll999l ISmail eHI\ h()Q2: 
ip 120021 K.K. Knudsen et al. in preparation). 
Best (2002) reported an excess number of submillimeter 
luminous galaxies (SMGs), detected at 850/im, in z ~ 
1 cluster fields and, based on lensing considerations, ar- 
gued that they were associated with the clusters rather 
than a background population. If the SMGs are indeed 
cluster members, they represent by far the most intense 
star bursts seen in clusters to date with rates of ^100 M 
yr _1 . Still, since none of the systems have been verified 
to lie at the cluster redshifts they could be background 
objects, perhaps lensed by the cluster potential. 

Lensing effects have been exploited by many groups 
to overcome the blank-field confusion limitations and 
detect faint SMGs (ffssoum < 2 mJy) at hig h red- 
shift iSmail. Ivison fc BlairJ Il997t iCowie et all 120021: 
iChapman et alJl2002fl : but these surveys have employed 
lower redshift clusters (z <; 0.5), since massive, relaxed 
clusters at low redshift are more efficient lenses than 
high-redshift clusters. Recent results from the Red - 
Sequence Cluster Survey (RCS; IGladders fc Yeell20f)^ . 
however, indicate that this is not always the case, and 
some higher redshift clusters appear to be excellent 

le nses. 

IGladders et a l. (2003) have found that the num- 
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ber of strong multiple optical arcs for z > 0.6 clus- 
ters seen in the RCS is larger than expected. This, 
along with the fact that the number of arcs in other 
medium- /high-redshift cluster surveys are several times 
larg er than that expected from standard ACDM cosmol- 
ogy iBartelmann et al.lll998t iZartiskv fc Gonzalezll2"003l: 
iDalal et all l2004|) . leads to the interpretation that at 
least at high redshift, there may be a subset of rich clus- 
ters which have m uch higher lensing cross- sections than 
predicted (but see iWambseanss et al.ll2004L for an alter- 
native explanation). These "super- lenses" could be due, 
for example, to an increase in cluster substructure with 
redshift, when clusters are observed in a younger and less 
relaxed state. 

To examine both the issues of super-lensing and and 
possible excess SMGs in high redshift clusters, we we 
have begun a survey at 850/^m of high redshift (0.6 
< z < 1.1) galaxy clusters from the RCS with SCUBA 
on the JCMT. The survey consists of two subsamples: 
five which show strong optical arcs and five control clus- 
ters of similar richness and redshift, which do not. Here 
we present the submillimcter data and follow-up imag- 
ing from a single super lensing cluster field RCSJ022434- 
0002.5 (hereafter RCS0224); the entire submm program, 
when complete, will be presented in a later paper (Webb 
et al., in preparation). This paper is organized as fol- 
lows. In §2 we outline the submillimcter and supple- 
mentary observations and data analysis; §3 presents the 
object catalogue; in §4 we discuss the properties of the 
counterparts to the submillimeter emission; §5 contains 
a discussion of the results. A flat, £l\=0.7 universe, with 
H =72 km/s/Mpc was assumed throughout. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. RCS0224 

RCS0224 inlanders. Yee. Ellin^sTirl [2002; 
Idadders et alJl2003(l was discovered within the RCS, a 
90-square degree optical imaging survey designed to find 
galaxy clusters up to z = 1.4. It is the most spectacular 
example of the RCS super lensing clusters, with two 
strong optical arcs originating from different background 
objects, one with z = 4.9, and a third arc was discov- 
ered in a deep HST image. Subsequent spectroscopic 
follow-up has confirmed the cluster redshift to be z = 
0.773. 

2.2. Submillimeter Data 

2.2.1. Observations and Reduction 

The cluster was simultaneously observed at 450/xm and 
850/im using the Su bmillimeter Common -User Bolome- 
ter Array (SCUBA) ijHolland et al.lll998fl on the James 
Clerk Maxwell Telescope (JCMT) over two observing 
runs in 2001 and 2002. We employed the jiggle-map 
mode of SCUBA which fills in the under-sampled sky by 
stepping the secondary mirror through a 64-point pat- 
tern in each nod position throughout each 128-sec in- 
tegration. With AC-coupled bolometers, chopping and 
nodding is still required in the submillimeter for proper 
sky removal and we chopped in RA with a 30" chop- 
throw. Such a small chop-throw has two advantages: it 
optimizes sky removal by sampling nearby sky, and the 
off-source chop remains on the array, creating a negative- 
positive-negative beam pattern that is used to improve 



the extraction of real sources. 

The map consists of two pointings, centrally offset from 
each other by 20" in RA and Dec and roughly centered on 
the cluster center, as determined from the optical data. 
The area is approxmately 6.3arcmin 2 though the depth 
varies with position in the map due to uneven coverage 
and differing bolometer sensitivities. A total integration 
time of 38.5ks was obtained on the center of the map, but 
the effective integration time is less at the edges where 
the sky is not evenly sampled by the array and where 
pointings do not overlap. 

The sky opacity was monitored through skydips ev- 
ery 1-2 hours and was stable throughout each individ- 
ual night, ranging from t 850 ~ 0.1-0.32 throughout the 
two runs. Pointing stability was also closely monitored: 
pointing observations on near-by point-sources were per- 
formed every 1.5 hours or after a major slew of the tele- 
scope. Calibration maps were acquired each night and 
calibration factors were found to be stable to within 10% 
of the standard gain values at 850/im. The calibration 
at 450/im is stable to ~ 30%. 

The data reduction was complicated by the presence 
of a strong correlated noise signal which occurred at the 
jiggle-map frequency (16 seconds) (C. Borys personal 
communication, 2004). This signal appears in the 2002 
data and affects between 1/3 and 2/3 of the bolometers 
at any given time. We reduced its effect in the follow- 
ing way. After the standard flat-fielding and extinction 
corrections, noise corrupted bolometers were identified 
through the presence of power in their fourier spectrum 
at the characteristic scale. Residual sky-flux was re- 
moved from all the bolometers by subtracting the me- 
dian sky level at each second, determined using only the 
non-corrupted bolometers. Because the noise on the re- 
maining corrupted bolometers is correlated it is possi- 
ble to reduce its effect through the subtraction of this 
correlated signal. This was done through simple mul- 
tiple linear regression techniques in which the expected 
noise signal in a given bolometer was predicted using all 
other corrupted bolometers and then removed. As a fi- 
nal step, noise spikes at > 3ct were iteratively removed 
from all bolometers and the bolometer time-streams were 
rcbinned to produce the final map. The removal of the 
correlated noise from the corrupted bolometers results 
in a overall reduction of the noise in the final map of 
~ 30%. A comparison of the non-corrupted 2001 data 
with the corrupted data from 2002, both corrected and 
non-corrected show that the presence of point sources 
in this map is robust, with no sources lost and no new 
sources introduced other than what is expected from an 
improved signal-to-noise ratio (S/N). We note, however, 
that this is not the case for shallower surveys where this 
noise signal can lead to significant chan ges in the number 
of re al and spurious sources detected l)Sawicki & Webbl 

Hoi. 

2.2.2. Source Detection 

Because of strongly varying noise properties (spatial 
and temporal), and the fact that one is generally work- 
ing at the detection limit of the data, deep imaging 
with SCUBA requires careful noise analysis, to assign 
significance to each detection, and to minimize spurious 
sources. To estimate the noise as a function of position 
on the image we employed two related techniques. First, 
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following th e metho d outlined in lEales ct al. ( 2000) and 
iWebb et all l|2003aD we generated 500 Monte Carlo sim- 
ulations of each bolometer time-stream, assuming Gaus- 
sian noise statistics, uncorrelated bolometers, and with- 
out introducing any signal. The reduction procedure was 
repeated on these simulated data, resulting in a set of 
final rebinned maps with no sources but with noise char- 
acteristics similar to the real data. A noise map is then 
simply the variance of the individual simulations. Sec- 
ond, using the real bolometer time-streams we employed 
a bootstrap algorithm to generate a variance map from 
500 realizations. The two techniques produced similar 
variance maps: the mean ratio between the two maps 
was ~ 1.0, and varied by <20% over the entire field. 
The significance levels quoted in Table^refer to the first 
method. 

Source extraction was performed u sing an iterative 
cleaning technique l|Eales et alJ I2000J) . An initial list 
of source positions was produced by convolving the 
raw map with the beam template, which includes the 
off-beams produced by the chop. This technique is 
simil ar to the templat e fitting used by other groups, 
(e.g.. iScott et alJ 12002'! and is advantageous over sim- 
ple smoothing since it incorporates information (position 
and flux) from the two negative off-sources associated 
with each real source. Beginning with an initial list of 
S/N > 3 detections, the map was iteratively cleaned of 
sources in 10% flux steps. Using these cleaned results, for 
each source all other sources were removed from the raw, 
unconvolved map, and the isolated source was again con- 
volved with the beam template. This provides improved 
flux densities, positions, and detection significances for 
each source since contamination from near-by confused 
sources is reduced (see Figure 1 and §3.1) 

The 450/im observations, which are taken simultane- 
ously by exploiting a dichroic beam splitter, suffer from 
a greatly increased sky opacity and a very unstable beam, 
and are therefore much more difficult to work with. In- 
stead of generating a separate source catalog from the 
450/im data we simply searched for 450/im emission from 
the previously detected 850/mr sources. We adopted a 
search radius of 10" around each 850/im position, allow- 
ing for a combination of the expected positional uncer- 
tainties at both wavelengths and we take as the flux mea- 
surement the nearest 450/im peak above 2.5ct within this 
radius. 

2.3. Supplementary Data 
2.3.1. Radio Imaging 

Approximately 5 hours of data were obtained on each 
of three days during December 2003 using the National 
Radio Astronomy Observatory's 8 (NRAO) Very Large 
Array in B configuration, exploiting a pseudo-continuum 
correlator mode to minimize bandwidth smearing. Data 
were recorded every 5 seconds in 3.25-MHz channels, 28 
in total, centered at 1.4 GHz, taking both left-circular 
and right-circular polarizations. 0137+331 was used for 
flux calibration and 0239+42 (0.72 Jy) was observed ev- 
ery hour for local phase/amplitude/bandpass calibration. 

Initial calibration and mapping were performed using 
standard AIPS tasks. The initial calibration was good 

8 NROA is operated by the Associated Universities Inc., under 
a cooperative agreement with the National Science Foundation. 



on most baselines, though we were able to improve the 
fidelity of the final image by running a series of self- 
calibration/imaging tasks, with a model comprising the 
central 820" x 820" field plus 11 non-contiguous fields 
containing bright sources. The final image has a series 
of north-south stripes near the brightest sources, typi- 
cal of equatorial fields, but in clean regions of the map 
the noise is around 15 /iJy beam -1 , where the beam is 
near-circular with fwhm ~ 4.5". 

2.3.2. Optical and Near-Infrared Imaging 

Multi-filter optical images of the field are available 
from various sources. The original discovery data from 
the RCS provide data in z' and R r using t he CFHT12K 
camera l|Gladders. Yee. fc Ellingsonl H)02) . Additional 
images in V and B have been obtained, also with the 
CFHT12K, as part of the four-color photometry follow- 
up of the CFHT part of the RCS survey. The details of 
the observatio ns and photomet r ic red uction techniques 
are provided in iGladders fc Yeel l)2005|) for the z' and R c 
data, and m lHsieh^t"ldTij2"005f) for the V and B images. 

J and if '-band images of the cen tral ~2'x2' were ob- 
tained using the PANIC camera l)Martini et all l2004|) 
on the Magellan Baade 6.5m telescope. PANIC is a 
1024x1024 pixel near-IR camera with a pixel scale of 
0.125" per pixel. Integrations totalling 54 min in K 1 were 
obtained on 21/12/2003 and 52 min in J on 24/10/2004, 
both under moderately good seeing and photometric con- 
ditions. Object finding and photometry in the K' image 
were performed using the program PPP (Yee 1991). The 
R c — K' colors were measured using identical angular 
aperture s izes based on the magnitude growth curve (see 
lYeelll99il for details). (J — K') colors, discussed in §5, 
were measured for the two brightest SMGs in a similar 
manner. 

3. SUBMILLIMETER RESULTS 
3.1. The Detected Sources 

The 850/im map is shown in Figure^^d the detected 
sources are listed in Table Q] We detect five SMGs at 
850/im above S/N > 3.0. Because the use of chopping to 
remove the sky level results in a map with a mean flux 
level of zero, the noise will be symmetric around zero 
and and the number of spurious sources on a SCUBA 
map may be estimated through the number of sources 
detected in an inverted map. After first removing all 
positive sources and their associated negative off-sources, 
we searched for negative sources at > 3.0 ct level on the 
inverted map. Only one negative source was found, at 
3. Oct near the northern edge of the map, indicating our 
source list is reliable at this level. 

To quantify the uncertainties in the submillimeter posi- 
tion and flux measurement we placed and recovered point 
sources of varying flux levels and positions into the im- 
age. Working with the real map will produce more real- 
istic results than working with the simulated data, since 
the effects of confusion are included through the pres- 
ence of the unresolved background. In order to avoid 
over-estimating the effects of confusion the five bright 
real sources (i.e., S/N > 3) were removed from the map 
before this analysis. Given our small field size and the 
high number of sources the addition of even a single ob- 
ject leads to an unrealistically large source density. On 
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the other hand, removing all the detected sources will 
lead to an underestimate of the positional uncertainties, 
since only confusion with sources below the confusion 
limit is considered. 

The 95th percentiles for the positional offsets for S/N 
of 3.5cr, 4.5cr, 5.5cr, and 6.5cr are 9.2", 8.4", 7.6", and 
7.6" respectively. The distributions in each bin are all 
centered at ~ 4.5" and skewed strongly to larger offsets, 
particularly in the lowest S/N bins. We note tha t these 
offsets are in good agreement with llvison et alJ l|2005f) 
who find a positional accuracy (la) of 1.5xFWHM/SNR 
for 5er sources, determined from the distribution of offsets 
from the objects to radio positions. The error on the flux 
measurement may be estimated from the signal-to-noise 
map, or through a comparison of the input flux to recov- 
ered flux in the simulations described above (more cor- 
rectly, from the output flux to the range in input fluxes 
which produce such a measured flux). It was found that 
these two methods compared well, indicating that the 
Monte Carlo simulations reasonably reproduce the noise 
statistics of the data. 

In Figure^b) we show the 450/im flux contours, over- 
laid on the 850/im map, and a correlation between the 
two wavelengths is apparent. Four of the five 850/im 
sources have usable 450/im coverage (the southern-most 
850/im source lies in the unreliable edge region of the 
450/im map). Three of these are detected at 450/um with 
> 3.0a, and the fourth with 2.5ct. Using the 450/im num- 
ber density on the map at these levels of significance we 
can estimate the probability that a 450/im detection is 
randomly aligned within 10" of an 850/im source. For 
the 450/im sources with S/N> 3 this probabilty is 8%, 
while the single S/N = 2.5cr has a probability of 16%. 
Given the actual measured offsets (which, barring one, 
are <J 6"), the probability of chance alignments for all 
four sources is < 6%. Thus, for these four sources we 
would expect less than one chance alignment and re- 
gard the 450/im detections as reliable counterparts to 
the 850/im emission. The detection of these four objects 
in both submillimeter bands is strong evidence that they 
are real and not spurious detections. 

In principle, because of the smaller beam size, the 
450/fin detection should provide a better measure of the 
true source position than at 850/im; in practice, the 
poorer quality of the 450/im data does not make this 
so. Simulations of the positional uncertainty at 450/im, 
similar to those discussed above for 850/im, indicate that 
at < 4a (the significance of these detections) the median 
positional uncertainty is roughly equal to that at 850/im 
but with higher dispersion and skewed heavily to large 
offsets (the 95% percentile is ~ 13"). 

3.2. Gravitational Lensing 

The magnifications were computed using a model that 
was derived to reproduce the strong lensing features ob- 
served for the cluster. The model consists of a smooth 
cluster component which was modeled as an isothermal 
ellipsoid with a core. The free model parameters are the 
central velocity dispersion, core radius, ellipticity, posi- 
tion angle and the position of smooth halo. In addition 
we included truncated isothermal spheres at the posi- 
tions of galaxies on the cluster Rc — z' color-magnitude 
relation. For these we assumed a Faber-Jackson relation 



and used the velocity dispersion of an L* galaxy as a free 
paramter in the model. 

We used deep optical HST observations 9 to model the 
lensing but found that the model cannot reproduce all 
strong lensing feature to great accuracy. The L* veloc- 
ity is not well constrained, but the best fit value of 180 
km/s is quite reasonable based on other studies (e.g., 
Hoekstra et al. 2000). The models prefer a small core 
radius. Despite the problems of reproducing all detailed 
lensing features, we found that the results for the mass, 
ellipticity and position angle of the smooth cluster com- 
ponent are robust. This is not surprising because these 
are well constrained by the high redshift arc system, for 
which we have confirmed the counter image spectroscop- 
ically. The orientation based on the strong le nsing model 
agree s well with that of the X-ray emission l|Hicks et alJ 
2004). The resulting mass estimate is also in good agree- 
ment with that determined from a weak lensing analysis 
of the HST image. 

Given the fact that most of the SMGs are found at 
relatively large distances from the cluster center, the de- 
rived strong lensing model is sufficiently accurate to de- 
rive their magnifications. In Table \I\we list the lensing 
magnifications these objects would experience assuming 
they lie behind the cluster. The magnifications are cal- 
culated assuming source plane redshifts of z — 1.5 - 4.0 
though in Table H we list only the two extremes. Four 
of the five of these objects lie at large angular distances 
from the cluster center and therefore suffer very small 
magnifications of ~ 1.1 for all plausible redshifts. Even 
source SMM-RCS0224.5 which lies close to the cluster 
center has a magnification of only 2.3 in the most ex- 
treme case. In conclusion, simple but accurate analysis 
shows that on average we do not expect the SMGs to be 
significantly magnified. 

4. THE RADIO AND NIR COUNTERPARTS AND 
REDSHIFT ESTIMATES 

To identify counterparts to the SMGs we employed the 
empi rical correlatio n between the radio and far-infrared 
flux l)Condonlll9 92^ observed in the local universe. Using 
the deep 1.4-GHz maps we searched for a radio detection 
within the submillimeter error radius (§3.1). Statistically 
this is also one of the least ambiguous methods of iden- 
tifying counterparts since the number density of radio 
sources on the sky is low enough that chance alignments 
between submillimeter and radio galaxies are almost neg- 
ligible. In Figures [3 and 01 we show the fC'-images sur- 
rounding each submm detection, with 1.4-GHz flux con- 
tours overlaid. Four of the five 850/^m sources have radio 
detections within the submillimeter search radius (§3.1) 
and all of these four have K' counterparts. The fifth 
source, SMM-RCS0224.4, does not have a secure coun- 
terpart at radio or optical/near- infrared wavelengths. 

The z'R c VB four-filter photometry data were used 
to generate photometric redshifts using an empirical 
training set technique, providing photometric redshifts 
with typical accuracy of ~ 10% for galaxies with R c <J 
23.0 (for details, see Hsieh et al. 2005). Two of the 

9 Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 5-26555. These 
observations are associated with program # 9135 
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TABLE 1 

Positions and Flux Densities of Submillimeter Sources 



Name 


R.A. 


Doc. 




850^tm 




450/xm 


offset: 850/im and 


fensing 




(J2000) 


(J2000) 


(mJy) 


S/N 


(mJy) 


S/N 


450/im pos.(") 


magniffcation a 


SMM-RCS0224.1 


02:24:34.29 


-00:03:28 


6.4 


7.0 


14.4 


3.5 


5.2 


1.10-1.17 


SMM-RCS0224.2 


02:24:33.56 


-00:03:58 


5.3 


4.6 


> 66 b 






1.06-1.10 


SMM-RCS0224.3 


02:24:29.79 


-00:03:01 


1.0 


3.3 


15.8 


3.0 


Cl-1 


1.05-1.09 


SMM-RCS0224.4 


02:24:32.69 


-00:03:46 


3.2 


3.1 


18.7 


2.5 


1.6 


1.06-1.11 


SMM-RCS0224.5 


02:24:33.00 


-00:03:35 


3.2 


4.1 


13.2 


3.0 


9.0 


1.61-2.59 



a For all sources but SMM-RCS0224.4 the lensing magnification was determined using the radio position (sec Tablc[2j. Two possible magnifications 
arc listed: the first assumes a source rcdshift of z — 1.5 and the second of z — 4.0. Note that this is simply the lensing calculated at the source 
position and none of these sources are confirmed as lensed objects. 

3<r upper-limit on the flux; the object lies at the edge of the usable 450/xm field. 



TABLE 2 
The Counterparts 



Name 


R.A. 


Dec 


1.4GHz flux 


offset from 




K' 


(Rc - K>) 


Zphot 




1.4GHz 


1.4GHz 


(MJy) 


85(tytm(") 






SMM-RCS0224.1 


02:24:33.815 


-00:03:33.00 


1652 ± 41 


8.7 


< o b 


16.92 ± 0.01 


4.63 ± 0.02 


0.66 ± 0.11 


SMM-RCS0224.2 


02:24:33.332 


-00:04:05.53 


68.4 ± 29 


8.3 d 


1.9 ± 1.0 


18.22 ± 0.04 


4.11 ± 0.05 


0.72 ± 0.08 


SMM-RCS0224.3 


02:24:29.972 


-00:03:04.10 


97.4 ± 29 


4.1 


0.85 ± 1.0 


20.81 ± 0.14 


~6.1 ± 0.7 




SMM-RCS0224.4 






< 48 e 




> 1.6° 








SMM-RCS0224.5 


02:24:32.867 


-00:02:31.34 


94.3 ± 28 


4.1 


0.65 ± 1.0 


19.38 ± 0.06 


4.6 ± 0.12 





a Using the empirical 1.4GHz-850/xm relation for high rcdshift submm-sclectcd galaxies from Chapman ct al. ( 2005). All estimates have uncertainties 
of Az ~ 1. 

^The substantial radio emission of this source, and the compact nature of the optical counterpart leads us to conclude that this object contains 
an AGN. As discussed in the text, this will result in an underestimate of the rcdshift and in this case a co mpletely non-physic al estimate of z < 0. 

c The photometric rcdshift is determined using the z'R c VB empirical training set technique presented in Hsich et al. (2005). Sources without a 
photometric rcdshift have not been detected in enough filters to enable an estimate. 

d The large positional offset could be due, in part, to blending with SMM-RCS0224.4 

e A 3(7 upper limit on the flux. 



four radio-detected objects (SMM-RCS0224.1 and SMM- 
RCS0224.2) have sufficiently accurate optical measure- 
ments to allow photometric redshift estimates. We were 
not able to obtain phtometric redshifts for the other 
two sources with radio counterparts. The sources SMM- 
RCS0224.3 is undetected in all our optical bands while 
SMM-RCS0224.5 is detected only marginally in the R c 
image. The redshift estimates are listed in Table along 
with the K' mag and Rc — K' colors. 

A detection at radio wavelengths has the added advan- 
tage that it provides a r edshift constraint, from the radio- 
subm m spectral index l|Dunne et al]l200"fl lYun fc Carillil 
2002). This method draws on the tight correlation be- 
tween FIR and radio flux for low redshift star form- 
ing galaxies and the steep slope of the submm side 
of the thermal dust spectrum (e.g., lEales et all 120001: 
llvison et jp Hig). There are, however, a number of 
substantial uncertainties in this technique including the 
assumed temperature of the dust and contamination of 
the radio flux by an AGN. Until recently the small num- 
ber of SMGs with spectroscopically determined redshifts 
made it difficult to verify the extension of this relation- 
ship to the high-redshift 850/im selected population. In 
Fig. U we show the 1.4-GHz-850/im spect ral index for the 
SMGs with spectroscopic redshifts from lGhanman et alJ 
( 2005) with two empirically determined relations from 
local starbursts overlaid. It is clear the SMG spectral in- 
dices are not well described by the empirical low-redshift 



relations at z ^ 1.4 and therefore we adopt the average 
relati on from the high redshift SMGs ijGhapman et alJ 
|2005|) to estimate redshifts for the four radio detected 
sources in the RCS0224 field (and a lower-limit for the 
fifth). These are listed in Table 2 and shown in Fig. 
We note that all have a substantial uncertainty of Az ~1. 

Source SMM-RCS0224.1 is particularly interesting. It 
is the brightest submm detection in the RCS0224 field 
and is identified with a bright and extended radio source 
with a large positional offset of 8.7". Since the positional 
error analysis outlined in §3.1 indicates that roughly 
5% of the detected sources should be recovered at >8" 
from their true position it is worrying that this is one 
of two radio identifications with >8" in this field. How- 
ever, given the extreme observed luminosity of the ra- 
dio source, it is statistically very unlikely that the radio 
source and submm so urce are not related (0.1% assuming 
the 1.4GHz counts of iBondi et al.l2003|) and we reiterate 
that the positional uncertainty could be underestimated 
(§3-1) 

The radio counterpart is asymetrically extended over 
7.0"x4.8", and centered on a highly compact, though re- 
solved, galaxy in the optical and NIR images. The 450/im 
emission (which should be resolved given the smaller 
450/im beam) shows hints of extension along the same 
position angle (see Figure ^) . The extreme radio flux 
of the source indicates the presence of an AGN: assum- 
ing the cluster redshift and a radio spectral index of 
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Fig. 1. — Left: The 850/im map of RCS0224, smoothed with a 14" Gaussian, with white indicating positive flux. Overlaid are the S/N 
contours (before cleaning) in 0.5 steps, beginning with 3. Oct and 1.0 steps above 4<r. The five sources listed in Table HI are marked. Sources 
2 and 4 are confused and separate into two sources during the 850/im cleaning process. This more easily apparent in the 450^tm map in 
which source 4 has been detected. Right: The greyscale 850/im map of RCS0224 with the 450/^m flux contours overlaid. The contours have 
been smoothed with an 8" Gaussian and are shown in steps of lmjy starting at 12mjy. Note that the 450^m FOV (denoted by the ellipse) 
is smaller than at 850/^m such that the outer ~20" of the 850/^m map does not have useful 450^m coverage, including the southern-most 
850/im. 



Socu^ - 7 yields a rest frame luminosity of 1.9xl0 24 W 
Hz -1 , well within the radio power regime dominated by 
AGN QCondon et all 12002]) . The asymmetric extended 
emission is suggestive of the morphology of a head-tail 
radio source, and specifically a Narrow Angle Tail galaxy 
(e.g.. ISiibrin fc de Bruvnlll998l) : a single, or parallel ra- 
dio jet, bent by bulk flows in the intercluster medium. 
The presence of a strong AGN can introduce uncertain- 
ties to the photometric redshift, however, this is unlikely 
to be significant in this case since the object is resolved 
optically. 

5. DISCUSSION 
5.1. SMG Source Counts 

A main goal of our submillimcter survey of RCS clus- 
ters is to determine whether some, or possibly all, z ~1 
galaxy cluster fields show an excess number of SMGs 
compared to that expected from the blank-field counts 
and simple lensing estimates. If so, this could be due ei- 
ther to enhanced gravitational lensing as explained in the 
§1 or to source count contamination by submillimctcr- 
bright cluster members. 

In Figure El we show the cumulative source counts to- 
ward RCS0224, with no correction made for lensing ef- 
fects, which we calculate to be negligible, but account- 
ing for the varying depths within the submillimeter map. 
The results from this single cluster field are consistent 
with the blank field counts and indicate an excess of 
sources due to either lensing or submillimeter-luminous 
cluster members only at the la level. Interestingly, the 
source counts are similar to thos e from other cluster sur- 
veys of iSmail etafl l)2002|) and iBestl l)2002j) which also 
lie systematically above the blank-fie ld counts at at least 
the la level (the exception to this is lGowie et al.l (2002) 
which is also a cluster program). The Smail et al. counts 
do not include two cluster members (cD galaxies) and 
have been corrected for lensing which can be substan- 
tial for low redshift clusters, whereas the higher redshift 
cluster fields of Best (and RCS0224) suffer a significantly 



lower average lensing magnification (§3.2). While the 
offset from the blank-field is small for each survey, taken 
together, these results suggest that source counts in clus- 
ter fields (at all redshifts) are enhanced, possibily due 
to cluster me mber contam i nation . This has been noted 
previously by iBorvs et alJ (|2003fl . The full RCS sub- 
millimeter survey will increase the statistics on cluster 
field counts, in particular at higher redshifts where the 
enhancement may be more pronounced, due to younger 
clusters and hierarchical processes (discussed below). 

5.2. Dusty Star Formation in RCS0224? 

Although the RCS0224 cumulate source counts are 
consistent with the blank-field measurements, this in it- 
self is not unambiguous evidence that all five SMGs are 
background systems; however, taken with our simple but 
adequate mass model (§3.2), it docs indicates that none 
of these systems are substantially lensed. Because of the 
scatter in the blank field source counts and the small 
area surveyed here, a scenario where some (but not all) 
of the SMGs are cluster members easily remains in-line 
with the expected source counts. 

The strongest evidence that some of the SMGs pre- 
sented here are cluster galaxies are the optical photomet- 
ric redshifts of two of the radio-detected SMGs, SMM- 
RCS0224.1 and SMM-RCS0224.2, which are in excellent 
agreement with the cluster redshift (Table[2J . The uncer- 
tainties on the optical photometric redshifts (e.g., Az ~ 
0.1) are substantially lower than for the radio-submm 
redshift estimate (Az ~1.0); for this reason whilst the 
radio-submm redshift estimates are discussed below we 
give them much less weight than the optical photometric 
redshifts for these two galaxies. 

The near-infrared counterparts of SMM-RCS0224.1 
and SMM-RCS0224.2 are shown in Fig. H and have 
colours and magnitudes (K' <; 18; see Table EJ sugges- 
tive of low-redshift objects. Their R c ~ K' colours are 
consistent with both the cluster red-sequence at this red- 
shift, (i? 702 -K) = 4.7 ijKodama et al.lll998|K as well as 
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Fig. 2. — K' images of the region surrounding each SMG. The images are centered on the SMG position with the approximate positional 
error indicated by a cross (8", ~95%). Overlaid are the 1.4GHz contours showing the radio detections. Source SMM-RCS0224.4 does not 
have a radio identification. The offset between the radio and optical positions for source SMM-RCS0224.2 is likely due to the low S/N of 
this object in the radio. East is to the left and North is up. 1.4GHz contours start at 40/^Jy and increase in 20/iJy steps up to 100/xJy, and 
then in 100/iJy steps. 



dusty star-bursts and radio gal axies at z ~ 1 l|Smail et alJ 
12004 iWaddineton et alll2000ft . The near-infrared color 
of SMM-RCS0224.1 is quite blue with (J-K)=l.l, simi- 
lar to other z ~ 1 radio galaxies (which this source cleary 
is with Sia-ghz = l-6mjy), and possibly indicating sub- 
stantial Ha and Nell contamination of the J-band flux. 
Source SMM-RCS0224.2 has (J - K) = 2.2, which if 
at z ~ 0.7 (the photometric rcdshift and cluster red- 
shift) sugges ts a dusty starburst, rather than an early- 
type galaxy l)Pozzetti fc Mannuccil l2000). 

The remaining two radio-detected SMGs are much 
fainter in the optical bands and thus do not have suffi- 
cient detections for phototmetric redshifts; however, the 



radio identification allows us to examine the R c — K' col- 
ors (our J imaging is too shallow for J — K). Source 
SMM-RCS0224.3 is very faint (K = 20.8) with an ex- 
tremely red color of at least (R c — K') ~ 6.1. This is 
consistent with a high-rcdshift z J> 1 galaxy and indeed 
with the counterparts of the z J> 2 blank-field SMGs 
which are often faint and red. Source SMM-RCS0224.5 
however, has a more moderate R c — K' color similar to 
SMM-RCS0224.1 and SMM-RCS0224.2. 

Redshifts estimated from radio-submm spectral indices 
have substantial uncertainties (Az ~ 1) and provide 
much poorer redshift constraints than the photometric 
redshifts. Nevertheless, since this is a frequently applied 
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Enlarged X-band images and radio contour maps of Fig|2]for SMM-RCS0224.3 (left) and SMM-RCS0224.5 (right), showing 



the faint X-band identifications and the close neighbors. 
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Fig. 4. — The radio-submm spectral index a = 0A2log(Ss50fj,m / Sia—GHz) as a functi on of redshift. The open diamonds correspond 
to submm-selected galaxies with 1.4GHz counterparts and spectros copic redsh ifts I Chapman et al. 2005 ). Also shown are two relations 
empirically determined using low-redshift starburst galaxies: solid line lYun & ; Carilli ( 2002) and dashed line lDunne et al] 120001) . The cluster 
redshift is denoted by the vertical solid line. The radio-submm redshift estimates for the four RCS0224 sources with radio detections are 
shown by the solid circles; these redshifts have been determined from the average relation of the submm-selected population (i.e., the open 
diamonds). The optical photometric redshift estimate of sources SMM-RCS0224.1 and SMM-RCS0224.2, arc shown by the solid stars and 
are connected to their corresponding radio estimate by dashed lines. We use the scatter of the lChapman et al J (2005) data to estimate the 
uncertainties of the radio-submm redshift redshift estimate which we plot as error bars for SMM-RCS0224.3 and SMM-RCS0224.5. For 
clarity, we have not plotted the errors for SMM-RCS0224.1 and SMM-RCS0224.2 which are similar to the other two sources. Note that the 
a value for SMM-RCS0224.1 implies a negative redshift which is likely due to the radio flux of this source being highly contaminated by 
emission from an AGN (see text). Source SMM-RCS0224.4 does not have a radio detection and the upper-limit on a, and the corresponding 
redshift limit are denoted by the errors. This limit carries the same uncertainty as the measured ratios. 



technique in the absence of spectroscopy, we discuss such 
estimates here. This method implies an unphysical red- 
shift of z < for SMM-RCS0224.1 which is almost cer- 
tainly due to the AGN contribution to the radio flux, 
in addition to the usual uncertainties. Two sources, 
SMM-RCS0224.3 and SMM-RCS0224.5 have measured 
ratios which are consistent with z < 1 objects. As dis- 
cussed above, however, the optical/NIR properties of 
SMM-RCS0224.3 do not support cluster membership but 
rather are similar to the background SMG population. 
Source SMM-RCS0224.5 on the other hand has a more 
moderate color and thus, while not ruling out a back- 
ground source, is in general agreement with the radio es- 
timate. The radio-to-submm redshift estimate of SMM- 
RCS0224.2 is z — 1.9, substantially higher than the op- 
tical/NIR estimate and while the radio-submm estimate 
is uncertain this raises the interesting possiblitity that 
the radio-submm galaxy and the optical counterpart may 
represent different galaxies, as there is a small offset (~ 
1 arcsec) between the radio and optical positions. 

In principle the 450/im-850/im flux ratios may also be 
used to isloate low-redshift SMGs but the large uncer- 
tainties on the 450/im flux measurement and the flatness 
of the observed ratio with redshift out to z ~ 2 (e.g. 



iWebb et al.l2003a|) makes this method of little use m our 
case. Three of the four 450/im detected SMGs are consis- 
tent with essentially all redshifts below z ~ 3, given a rea- 
sonable thermal spectral energy distribution (SED). The 
fourth source, the strong radio source SMM-RCS0224.1, 
has a ratio which appears to place it at z > 3; however 
in this case it is likely the 850/zm flux is also boosted 
by the AGN, though to a lesser extent than the radio 
flux, and therefore the 450/im-850/im flux ratio cannot 
be described by a simple thermal SED. 

To summarize, sources SMM-RCS0224.1, SMM- 
RCS0224.2 are likely cluster members, based on their 
optical photometric redshifts, K magnitudes and opti- 
cal/NIR colors. Sources SMM-RCS0224.3 and SMM- 
RCS0224.4 are more consistent with belonging to the 
background population of SMGs, though they are not 
highly lensed by the cluster. Very little constraint can 
be placed on the redshift of SMM-RCS0224.5; there are 
not enough data for an optical photometric redshift es- 
timate and while its radio-submillimeter ratio suggests 
a z < 1 galaxy, there is substantial uncertainty in this 
technique, as evidenced by Fig.0] Spectroscopic redshift 
measurements are clearly required to resolve this issue, 
as well as improved positional information for the SMGs 
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Fig. 5. — The cumulative source counts for RCS0224 (large solid stars), assuming none of the detected objects have been lensed by the 
cluster. Also shown are the source counts determined from blank field surveys and cluster lensing surveys, corrected for lensing (Smail et 
al. and Cowie et al.) and not corrected (Best). Note that 3/4 of the cluster surveys (Smail, Best, this work) lie systematically above the 
blank field counts at the la level. As mentioned in the text, this may be an indication of unrecognized submillimeter-luminous cluster 
members enhancing the counts, even at low redshift. 



with large offsets to their radio identifications. The lat- 
ter is particularly important to rule out the possibility of 
single galaxy, or cluster-substructure, lensing. 

If two of the SMGs in the RCS0224 field are indeed 
dusty starbursts and/or AGN within the cluster this has 
interesting implications for the evolution of galaxy clus- 
ters. Thus far, star formation rates and dust levels great 
enough to produce strong enough submillimeter emis- 
sion to be seen by SCUBA (~100 M Q yr _1 , but very 
dependent on temperature) have not been unambigously 
detected in clusters. This could, in part, be a selection 
effect: in addition to the fact that very few high redshift 
galaxy clusters have been searched for star formation, 
such objects may be so obscured that they go unrecog- 
nized in optical studies, as was the case in the field un- 
til recently. Submillimeter surveys are sensitive to much 
colder dust than mid-IR observations with ISO, and thus 
may detect fundamentally different objects. In fact, a 
deep SCUBA image has been obtained for A1689 but 
none of the ISO sources are coincident with the 850/im- 
selected objects (Knudsen et al. in preparation). 

We can estimate the star formation rate such systems 



would have if they lie at z = 0.773 through an extrap- 
olation from Ssso/im to Lpm and using the relation of 
iKennicufrq l)1998f) . We ad opt a temperature of 40K and 
dust emissivity j3 = 1.3 l|Dunne et al.1 120 00). Ignoring 
SMM-RCS0224.1, which is likely AGN-dominated, we 
find a range in star formation rates of ~ 250-430 M Q 
yr _1 , for our range in 850/zm flux levels. We note, how- 
ever, that this estimate is h ighly dependent on temper- 
ature fseej Eales et alJ l2000. for discussion) and varying 
the dust temperature from 30-50K (a relatively modest 
range) results in SFR estimates from 90-900Mq yr . 

We note that both of the clusters in the literature with 
subst antial 15/im emission (jFadda et alJl2000HDuc et alJ 
2002) show signs of recent or ongoing cluster mergers, 
as do many of the clusters with excess radi o sources at 
high redshift and within the local universe jSmail et alJ 
ll999tlMiller fc Owenl2003llGiacintucci et al.l2004fl . Cur- 
rently there is no evidence that RCS0224 is undergoing 
a merger or large accretion, other th an a shallow X-ra y 
map showing hints of substructure (Hi cks et al.l 12004). 
but this is simply a reflection of the limited amount of 
data acquired for this cluster so far. Interestingly, an 
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increase in cluster merging at z ~ 1 also provides an 
explanation for the super -lensing phenomen a, of which 
RCS0224 is an example. iTorri et al] (|2004|) show that 
during a merger the lensing cross-section of a cluster is 
temporarily enhanced to such a degree that this process 
can account for the excess multiple arcs seen in optical 
surveys. In such a picture, the presence of strong optical 
arcs may mark merging clusters in which we may also 
find enhanced star formation and AGN activity. 

6. CONCLUSIONS 

We are undertaking an 850/im imaging survey of z~l 
galaxy clusters selected from the Red-Sequence Cluster 
Survey. Here we present the submillimeter data on the 
first completed cluster, RCS0224, and follow-up imag- 
ing for counterpart determination. We detect five SMGs 
within the field which our mass model shows are not 
strongly lensed, and measure cumulative number counts 
that are in general agreement with that expected from 
other blank-field surveys. Interstingly, this field and two 
of the three other cluster surveys lie systematically above 
the blank-field counts, all at roughly la, and taken to- 
gether suggest that SMG counts in cluster fields may be 
enhanced by unrecognized submillimeter luminous clus- 
ter members. A larger sample of clusters will provide 



tighter statistics. 

All five SMGs in the RCS0224 field are detected at 
1.4GHz or 450^m, or both, and we have located NIR 
counterparts for all the radio detections. Based on 
accurate optical photometric redshifts, and magnitude 
and color considerations, we find that two of the five 
SMGs (SMM-RCS0224.1, SMM-RCS0224.2) are consis- 
tent with being cluster member while another two ob- 
jects are more likely background objects. The redshift of 
the last source (SMM-RCS0224.5) is poorly constrained; 
while its radio-submillimctcr flux ratio is consistent with 
the cluster redshift the large uncertainty in this tech- 
nique makes this measurement unreliable. Further data, 
such as spectroscopic redshifts and improved submillime- 
ter positions, will verify this picture. 



We thank the anonymous referee for a careful read- 
ing of the manuscript. Research by Tracy Webb is sup- 
ported by the NOVA Postdoctoral Fellowship program. 
The JCMT is operated by the Joint Astronomy Cen- 
ter on behalf of the Particle Physics and Astronomy Re- 
search Council in the UK, the National Research Council 
of Canada, and the Netherlands organization for Scien- 
tific Research. 



REFERENCES 



Bartelmann, M., Huss, A., colberg, J.M., Jenkins, A., & Pearce, 

F.R., 1998, A&A, 330, 1 
Best, P.N., van Dokkum, P.G., Franx, M., & Rottgering, H.J. A. 

2002, MNRAS, 330, 17 
Best, P. 2002, MNRAS, 336, 1293 

Blain, A., Chapman, S.C., Smail I., & Ivison, R., 2004, ApJ, 611, 
52 

Bondi, M., et al., 2003, A&A, 403, 857 

Borys, O, et al. 2003, MNRAS, 344, 385 

Butcher, H., & Oemer, A., 1984, ApJ, 285, 426 

Chapman, S.C., Scott, D., Borys, O, & Fahlman, G.G. 2002, 

MNRAS, 330, 92 
Chapman, S.C., Blain, A.W., Smail, I., & Ivison, R., 2005, Apjin 

press, astro-ph/0412574 
Coia, D., et al., '2WA, A&A, in press, astro-ph/0310317 
Condon, J. J., 1992, ARA&A, 30 575 

Condon, J. J., Cotton, W.D, & Broderick, J.J., 2002, AJ, 124, 675 
Cowie, L.L., Barger, A.J, & Kneib, J.-P, AJ, 123, 2197 
Dalai, N., Holder, C, & Hennawi, J.F., 2004, ApJ, 609, 50 
Due, P.-A., et al., 2002, A&A, 382, 60 

Dunne, L., Clements, D.L., & Eales, S.A., 2000, MNRAS, 319, 813 

Dwarakanath, K.S. & Owen, F.N. 1999, ApJ, 118, 625 

Eales, S., Lilly, S., Webb, T., Dunne, L., Gear, W., Clements, D., 

& Yun, M., 2000, AJ, 120, 2244 
Edge, A.C., Ivison, R.J., Smail, I., Blain, A.W., & Kneib, J.-P. 

1999, MNRAS, 306, 599 
Fadda, D., Elbaz, Due, P.-A., Flores, H., Franceschini, A., 

Cesarsky, C.J., & Moorwood, A.F.M, 2000, A&A, 361, 827 
Giacintucci, S., Venturi, T., Bardelli, S., Dallacasa, D., & Zucca, 

E., 2004, A&A, 419, 71 
Gladders, M.D., Yee, H.K.C., & Ellingson, E., 2002, AJ, 123.1 
Gladders, M.D., Hoekstra, H., Yee, H.K.C, Hall, P.B., & 

Barrientos, L.F., 2003, ApJ593, 48 
Gladders, M.D. & Yee, H.K.C. 2005, ApJS, 157,1 
Hicks, A., Ellingson, E., Bautz, M., Yee, H.K.C, Gladders, M., 

& Garmire, G., 2004, Advances in Space Research: Clusters 

of Galaxies: New Insights from XMM-Newton, Chandra and 

INTEGRAL, in press, |astro-ph/0410167| 
Holland, W.S., Cunningham, U.K., Gear, W.K., Jenness, T., 

Laidlaw, K., Lightfoot, J.F., & Robson, E.I., 1998, SPIE, 3357, 

305 

Hoekstra, H, Franx, M., & Kuijken, K., 2000, ApJ, 532, 88 



Hsieh, B.C., Yee, H.K.C, Lin, H., & Gladders, M.D., 2005, ApJS, 

in press, astro-ph/0502157 
Ivison, R.J. et al. 2002, MINKAS, 337, 1 

Ivison, R.J., 2005, in The Spectral Distributions of Gas 
Rich Galaxies: Confronting Models with Data, in press, 
astro-ph/0412084 

Kennicutt, K., 1998, Ap.J, 198, 511 

Kodama, T., Arimoto, N., Barger, A. J., & Aragoh-Salamanca, A. 

1998, A&A, 334, 99 
Martini, PI, Persson, S.E., Murphy, D.C, Birk, C, Shectman, S.A., 

Gunnels, S.M., & Kock, E., 2004, Proc. SPIE: Astronomical 

Telescopes and Instrumentation, in press, astro-ph/0406666 
Miller, N.A., & Owen, F.Z., 2003, AJ, 125, 2427 
Poggianti, B.M., Smail, I., Dressier, A., Couch, W.J., Barger, A. J., 

Butcher, H., Ellis, R.S. & Oemler, A. 1999, ApJ, 518, 576 
Pozzetti, L & Mannucci, F, 2000, MNRAS, 317, 17 
Sawicki, M. & Webb, T.M.A., 2004, ApJ, in press 
Scott, S. et al. 2002, MNRAS, 331, 817 
Sijbring, D. & de Bruyn, A.G., 1998, A%A, 331, 901 
Smail, I., Ivison, R.J., & Blain, A.W., 1997, ApJ, 490, 5 
Smail, I., Morrison, C, Gray, M.E., Owen, F.N., Ivison, R.J., 

Kneib, J.-P., & Ellis R.S. 1999, ApJ, 525, 609 
Smail, I., Ivison, R.J., Blain, A.W., & Kneib, J.-P., MNRAS, 331, 

495 

Smail, I., Chapman, S.C, Blain, A.W., & Ivison, R.J., ApJ, 616, 
71 

Torri, E., Meneghetti, M., Bartelmann, M., Moscardini, L., Rasia, 

E., & Tormen, C, 2004, MNRAS, 349, 476 
van Dokkum, P.G., Franx, M., Fabricant, D., Kelson, D.D., & 

Illingworth, CD. 1999, ApJ, 520, L96 
van Dokkum, P.C, Stanford, S.A., Holden, B.P., Eisenhardt, P.R., 

Dickinson, M., & Elston, R. 2001, ApJ, 552, L101 
Waddington, I., Windhorst, R.A., Dunlop, J.S., Koo, D.C, & 

Peacock, J. A., 2000, MNRAS, 317, 801 
Wambsganss, J., Bode, P., & Ostriker, J. PI, 2004, ApJ, 606, L93 
Webb, T.M.A., Eales, S.A., Lilly, S.J., Clements, D.L., Dunne, L., 

Gear, W.K., Ivison, R.J., Flores, H., & Yun, M., 2003, 587, 41 
Yee, H.K.C, 1991, PASP, 103, 396 
Yun, M. & Carilli, C, 2002, ApJ, 568, 88 
Zartisky, D. & Gonzalez, A.H., 2003, ApJ, 584, 691 



